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ABSTRACT 

We present detailed NLTE Ti I calculations in model atmospheres of cool 
dwarf and giant stars. A fully self-consistent NLTE treatment for a Ti I model 
atom with 395 levels and 5279 primary bound-bound transitions is included, and 
we discuss the implication of departures from LTE in this atom for the strengths 
of Ti I lines and TiO molecular bands in cool star spectra. We show that in 
the atmospheric parameter range investigated, LTE is a poor approximation to 
Ti I line formation, as expected from the low collisional rates in cool stars. The 
secondary effects of Ti I overionization on the TiO number density and the TiO 
molecular opacities, however, are found to be negligible in the molecular line 
forming region for the relatively small parameter range studied in this paper. 

1. Introduction 

The opacity in the outer layers of the atmospheres of M stars are dominated by a 
small number of very strong molecular compounds (H 2 0, TiO, H 2 , CO, VO). Most of the 
hydrogen is locked in molecular H2, most of the carbon in CO; and H2O, TiO, and VO 
opacities define a pseudo-continuum covering the entire flux distribution of these stars. The 
optical "continuum" is due to TiO vibrational bands which are often used as temperature 
indicators for these stars. The indicators may be the depth of the bands relative to the 
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peaks between them; or the depth of the VO bands; or the strength of the atomic lines 
relative to the local "continuum" ; or even the strength of the infrared water bands; all of 
these depend on the strength of the TiO bands and the amount of flux-redistribution to 
longer wavelengths caused by the blocking of the flux in the optical and near-IR by TiO. 
Departures from LTE of the Ti I atom, and thus indirect changes in the concentration of the 
important TiO molecule, could therefore have severe and measurable consequences on the 
atmospheric structure and spectra of cool stars (Johnson, 1994; Gustafsson & Jorgensen, 
1994). Departures from LTE in the photospheres of cool stars have been investigated 
in detail in the line transfer and dissociative equilibrium of H 2 in red giant atmospheres 
(Lambert & Pagel, 1968; Vernazza, Avrett, & Loeser, 1981; Anderson, 1989). NLTE effects 
of CO have been investigated in the Sun (Thompson, 1973; Ayres & Wiedemann, 1989) and 
in red giant atmospheres (Carbon, Milkey, & Haesley, 1976; Wiedemann & Ayres, 1991; 
Wiedemann et al., 1994). These studies led to the conclusion that even after the passage of 
pulsationally driven shocks the recombination of H2 and CO must proceed rapidly, leading 
to number densities very close to LTE values. Pavlenko et al. (1995) and Carlsson et al. 
(1995) showed that the equivalent widths of the Li I resonance lines are affected by NLTE 
effects and have derived NLTE curves of growth for a range of model parameters. Andretta, 
Doyle, & Byrne (1997) have performed NLTE calculations for Na I using a method similar 
to that of Carlsson et al. (1995). NLTE effects on the excitation potential and spectral 
lines of the much larger Ti atom have been analyzed in late-type giants by Ruland et al. 
(1980); Brown, Tomkin, & Lambert (1983). 

Due to their very low electron temperatures, the electron density is extremely low in M 
stars; the absolute electron densities are even lower than found in low density atmospheres, 
such as those of novae and SNe. Collisions with particles other than electrons, e.g., H 2 or 
helium, are not as effective as electron collisions in restoring LTE, both because of their 
smaller cross-sections and their much smaller thermal velocities. Therefore, collisional rates 
which tend to restore LTE, could be very small in cool stars. This in turn could significantly 
increase the importance of NLTE effects in M stars when compared to, e.g., solar type stars 
with much higher electron densities and temperatures. 

The NLTE treatment of molecules such as H2O and TiO which have several million 
transitions is a formidable problem which requires an efficient method for the numerical 
solution of the multi-level NLTE radiative transfer problem. Classical techniques, such as 
the complete linearization or the Equivalent Two Level Atom method, are computationally 
prohibitive for large model atoms and molecules. Currently, the operator splitting or 
approximate A-operator iteration (ALI) method (e.g., Cannon, 1973; Rybicki, 1972, 1984; 
Scharmer, 1984) seems to be the most effective way of treating complex NLTE radiative 
transfer and rate equation problems. Variants of the ALI method have been developed to 
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handle complex model atoms, e.g., Anderson's multi-group scheme (Anderson, 1987, 1989) 
or extensions of the opacity distribution function method (Hubeny & Lanz, 1995). However, 
these methods have problems if line overlaps are complex or if the line opacity changes 
rapidly with optical depth, a situation which occurs in cool stellar atmospheres. The ALI 
rate operator formalism, on the other hand, has been used successfully to treat very large 
model atoms such as Fe directly and efficiently (cf. Hauschildt & Baron, 1995; Hauschildt 
et al., 1996; Baron et al, 1996). 

In this paper we discuss NLTE effects of Ti I in fully self-consistent models for a few 
representative M dwarf and M giant model atmospheres and spectra. 

2. Methods and Models 

In order to investigate the importance of Ti I NLTE effects on the formation of 
cool star spectra, a full NLTE model calculation is required. This means that the 
multi-level NLTE rate equations must be solved self-consistently and simultaneously with 
the radiative transfer and energy equations, including the effects of line blanketing and of 
the molecular equation of state. For the purpose of this analysis we use our multi-purpose 
stellar atmosphere code PHOENIX originally developed for the modeling of the expanding 
atmospheres of novae and supernovae, and adapted to conditions found in cool stars by 
Allard and Hauschildt (1995, hereafter AH95). PHOENIX (version 8.0, Hauschildt et al., 
1996; Hauschildt, Baron, & Allard, 1997) uses a spherical radiative transfer for giant models 
(log(g) < 3.5), and an equation of state (EOS) including more than 600 molecules. In 
all models TiO + (and ZrO + ) species are included in the equations of state according to 
partition functions published by Gurvich & Glushko (1982) for a consistent ionization 
equilibrium of Ti. For the strongest ca. 4 x 10 6 atomic & molecular lines, we use detailed 
depth-dependent Voigt profiles with improved damping constant computation (Schweitzer, 
1995; Schweitzer et al., 1996), and Gaussian profiles for an additional 10 7 much weaker 
lines. In addition, we include ca. 2000 photo-ionization cross sections for atoms and ions 
(Verner & Yakovlev 1995). These and related improvements, as well as the resulting LTE 
"NextGen" model grids, are described in detail by Schweitzer (1995); Schweitzer et al. 
(1996) and Allard & Hauschildt (in preparation). 

Both the NLTE and LTE lines for atoms and molecules are treated with a direct 
line-by-line method. We do not use pre-computed opacity sampling tables, and thereby allow 
for computation of pressure-dependent Voigt profiles, changes in the isotope ratios, and 
avoid possible interpolation effects. We dynamically select the relevant lines from master 
line lists at the beginning of each iteration and sum the contribution of every line within a 
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search window to compute the total line opacity at arbitrary wavelength points. The latter 
is crucial in NLTE calculations in which the wavelength grid is both irregular and variable 
(from iteration to iteration due to changes in the physical conditions). This approach 
also allows detailed and depth dependent line profiles to be used during the iterations. 
To make this method computationally efficient, we employ modern numerical techniques, 
e.g., vectorized and parallel block algorithms with high data locality (Hauschildt, Baron, 
& Allard, 1997), and we use on high-end workstations or supercomputers for the model 
calculations. In the calculations we present in this paper, we have set the micro-turbulent 
velocity £ to 2kms" 1 . We include LTE lines (i.e., lines of species that are not treated in 
NLTE) if they are stronger than a threshold r = Xi/ K c — 10~ 4 , where xi is the extinction 
coefficient of the line at the line center and k c is the local b-f absorption coefficient. 

PHOENIX is a full multi- level NLTE code, i.e., NLTE effects are considered self- 
consistently in the model calculations, including the temperature corrections. The 
temperature corrections and the convection are treated as described in Allard & Hauschildt 
(1995). Hauschildt & Baron (1995) have extended the numerical method developed by 
Hauschildt (1993) for NLTE calculations with a very detailed model atom of Fe II. In this 
section we describe how we apply this technique to a detailed Ti I model atom from the list 
of NLTE species already available in the model calculations. 



2.1. NLTE Calculational Method 

The large number of transitions of the Ti I ion that have to be included in realistic 
models of the Ti I NLTE line formation require an efficient method for the numerical 
solution of the multi-level NLTE radiative transfer problem. As already mentioned, the 
Ti I model atom described here includes more than 5200 individual NLTE lines plus a 
large number of weak background transitions. Classical techniques, such as the complete 
linearization or the Equivalent Two Level Atom methods, are computationally prohibitive. 
In addition, we are also modeling moving media (e.g., stellar winds, novae and supernovae), 
therefore, approaches such as Anderson's multi-group scheme or extensions of the opacity 
distribution function method (Hubeny & Lanz (1995)) cannot be applied. Again, simple 
approximations such as the Sobolev method, are very inaccurate in problems in which lines 
overlap strongly and make a significant continuum contribution (important for weak lines), 
as is the case for nova (and SN) atmospheres (Hauschildt et al., 1996; Baron et al., 1996, 
cf.). 



We use, therefore, the multi-level operator splitting method described by Hauschildt 
(1993). This method solves the non-grey, spherically symmetric, special relativistic equation 
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of radiative transfer in the co-moving (Lagrangian) frame using the operator splitting 
method described in Hauschildt (1992). Details of the method are also described in 
Hauschildt & Baron (1995), so we give here only a summary of the method. 

Even with highly effective numerical techniques, the treatment of possibly more than 
one million NLTE lines poses a significant computational problem, in particular in terms of 
memory usage. In addition, most lines are very weak and do not contribute significantly 
to the radiative rates. However, together, they can influence the radiation field from 
overlapping stronger transitions and should be included as background opacity. Therefore, 
we separate the stronger "primary" lines from the weaker "secondary" lines by defining a 
threshold in log(gf ), which can be arbitrarily changed. Lines with g /-values larger than 
the threshold are treated in detail, i.e., they are fully included as individual transitions in 
the radiative transfer (assuming complete redistribution) and rate equations. In addition, 
we include special wavelength points within the profile of the strong primary lines. 

The secondary transitions are included as background NLTE opacity sources but 
are not explicitly included in the rate equations. Their cumulative effect on the rates 
is included, since the secondary lines are treated by line-by-line opacity sampling in the 
solution of the radiative transfer equation. Note that the distinction between primary and 
secondary transitions is only a matter of convenience and technical feasibility. It is not a 
restriction of our method or the computer code but can be easily changed by altering the 
appropriate input files. As more powerful computers become available, all transitions can 
be handled as primary lines by simply changing the input files accordingly. We do not pose 
additional thresholds such as the energy or the statistical weight of the lower level of a 
line. However, we include in the selection process only observed lines between known levels 
in order to include only lines with well known gf- values. All predicted lines of Kurucz 
are included as secondary lines (see Hauschildt & Baron, 1995) for completeness. For all 
primary lines the radiative rates and the "approximate rate operators" (Hauschildt, 1993) 
are computed and included in the iteration process. 



2.2. The Ti I model atom 

To construct the Ti I model atom we have selected the first 34 terms of Ti I. We 
include all observed levels that have observed b-b transitions with log (gf) > —3.0 as NLTE 
levels where g is the statistical weight of the lower level and / is the oscillator strength of 
the transition. This leads to a model atom with 395 levels and 5279 primary transitions 
treated in detailed NLTE. That is, we solve the complete b-f & b-b radiative transfer 
and rate equations for all these levels including all radiative rates of the primary lines. A 
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Grotrian diagram of this model atom is shown in Fig. 1. In addition, we treat the opacity 
and emissivity for the remaining 0.8 million weak secondary b-b transitions in NLTE, if 
one level of a secondary transition is included in the model. A detailed description of the 
numerical method is presented in Hauschildt & Baron (1995). 

Photo-ionization and collisional rates for Ti I are not yet available. Thus, we have 
taken the results of the Hartree Slater central field calculations of Reilman & Manson 
(1979) to scale the ground state photo-ionization rate and have then used a hydrogenic 
approximation for the energy variation of the cross-section. Although they are only very 
rough approximations, the exact values of the b-f cross- sect ions are not important for 
the opacities themselves which are dominated by known b-b transitions of Ti I and other 
species. They do, however, have an influence on the actual b-f rates but this remains 
unimportant for the computational method used in this work. 

While collisional rates are important in hotter stellar atmospheres with high electron 
densities, they remain nearly negligible when compared to the radiative rates for the low 
electron densities found in cool stars. We have approximated bound-free collisional rates 
using the semi-empirical formula of Drawin (1961). The bound-bound collisional rates are 
approximated by the semi-empirical formula of Allen (1973), while the Van Regemorter's 
formula (1962) was used for permitted transitions. 

A more accurate treatment of this model atom requires the availability of more 
accurate collisional and photo-ionization rates for Ti I. In the present calculations we 
have neglected collisions with particles other than electrons because the cross-sections are 
basically unknown. Additional collisional processes would tend to restore LTE, therefore, 
the NLTE effects that we obtain in our calculations should be maximized. 



3. Results 

We have computed a small number of both giant and dwarf models to investigate the 
effects of Ti NLTE on the structure and the spectra of cool stars. The giant models are 
computed using the spherically symmetric, static mode of PHOENIX, while the dwarf models 
use the plane-parallel mode. For the three giants we use a gravity of log(g) = 1.0, solar 
abundances and the following T c g = 3200 K, 3600 K, and 4000 K. All giant models have 
radial extensions of less than 10%. In addition, we computed two NLTE dwarf models with 
log((?) = 5.0, solar abundances and T c g = 2700 K and 4000 K for comparison. All models 
include the Ti I NLTE treatment as discussed above as well as the standard PHOENIX 
equation of state (NLTE mode) and additional LTE background lines (about 10-15 million 
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atomic and molecular lines). The NLTE effects are included in both the temperature 
iterations (so that the structure of the models includes NLTE effects) and all radiative 
transfer calculations. With the exception of the hottest models, the majority of the line 
opacity comes from the LTE TiO and water wapor lines in most layers of the models, 
however, the Ti I lines are an important source of opacity in some of the outer laters of the 
models, in particular in the hotter models. For Ti I primary lines we use 5 to 1 1 wavelength 
points within their profiles. This procedure typically leads to about 150,000 wavelength 
points for both the model iteration and the synthetic spectrum calculations. 

3.1. Departure coefficients for Ti I 

In Figs. 2 — 4 we show the departure coefficients of Ti I for the three giant models. 
Figures 5 and 6 show the corresponding results for the two dwarf models. The optical depth 
scale, r st d, is measured in the b-f and f-f continuum at 1.2/im. Note that the atmospheres 
are extremely non-grey, some Ti I lines form at r st d ~ 10~ 6 . Therefore, we show in the 
following graphs a large dynamical range of standard optical depth. 

The departure coefficients in the dwarf models are significantly smaller than in the 
giant models, by several orders of magnitude. This is mostly an ionization effect, in the 
dwarf models Ti I and Ti molecules (mainly TiO) have higher concentrations relative to 
Ti II than in the giant models (see below). The electron temperatures in the outer regions of 
the dwarfs and giants are comparable. Therefore, the departure coefficients bi are distinctly 
different since their definition explicitly includes the concentrations of both electrons and 
Ti II (cf. Mihalas, 1978, p. 219). 

For the giant models, the T ef[ = 3600 K model has the smallest spread of the departure 
coefficients compared to the 3200 K and 4000 K models, indicating the smallest "internal" 
NLTE effects for this model. The range of the b^s in all giant models is about a factor of 10 
at T s td — 10~ 6 , but the line forming regions are farther inward, roughly between r st d = 10~ 4 
and T st d = 1- For most of the Ti I levels, the departure coefficients are less than unity, 
however at lower effective temperatures more levels show bi > 1 in the outer part of the 
atmosphere. Typically, b; L < 1 can be associated with an overionization of Ti I relative to 
the LTE state. However, this is valid only if the electron densities of the NLTE and LTE 
cases are very similar, otherwise the behavior is more complicated (Hauschildt et al., 
1996). In the models presented here, the electron densities of the LTE and NLTE models 
are basically identical because we have intentionally treated the important electron donors 
(mostly alkali and earth-alkali metals) in LTE. 
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The spread of the departure coefficients is much larger in the dwarf models, in 
particular for the cooler model. However, the large electron pressures in the dwarfs will 
compensate for the effect of the departures from LTE and result in smaller changes in the 
Ti I line profiles, as we discuss below. 

3.2. NLTE effects on balance of Ti I, II and Ti-molecules 

To analyze the effects of Ti 1 NLTE on the formation of the important TiO molecule, the 
major opacity source in the optical spectra of cool stars, we plot the relative concentration 
Pi/Pgas °f a variety of Ti ions and molecules for both the LTE and the NLTE cases in 
Figs. 7 — 11. The LTE figures were constructed by using the NLTE structure of the model 
atmosphere but setting all hi = 1. The LTE and NLTE structures are very similar out 
f° Tstd ~ 10~ 7 , in the very outer regions the NLTE models are up to 300K cooler than 
their LTE counterparts due to increased line cooling. Therefore, we have used the NLTE 
structures exclusively in these figures in order to avoid the effects of the cooler outer layers 
on the plots. 

The main differences between the giant and the dwarf models are caused by higher 
pressures in the dwarf models compared to the giants. This results in Ti II being more 
important in the giants than in the dwarfs (compare Figs. 9 and 11). The higher pressures 
in dwarfs also cause much higher concentrations of Ti-molecules in dwarfs than in giants of 
comparable effective temperatures. 

Comparison of the LTE and NLTE results for the same effective temperatures shows 
that for the giants the effects on the concentration of Ti species is relatively small for 
T cS = 3200 K. Only for r std < 10" 4 do NLTE effects on the EOS become noticeable. The 
main effect is that the concentration of Ti I is smaller, whereas the concentration of Ti II 
and TiO + are slightly higher than in the LTE model. The concentration of TiO is hardly 
changed by NLTE effects in this particular model, only in the outermost optically very thin 
regions does the TiO concentration drop, with a maximum change of about a factor of 2. 

The results are similar for the T e g = 3600 K giant model. Here, NLTE effects prevent 
recombination of Ti II to Ti I in the very outer parts of the atmosphere. However, the 
optical depths there are very small and the effects on the spectrum due to the EOS NLTE 
effects are quite small. As in the T e s = 3200 K model, the TiO concentration is reduced 
but the TiO + concentration is increased in the regions where NLTE effects are important 
for the EOS. The T cS = 4000 K giant model reacts similarly to the 3600 K model, but 
at this effective temperature Ti II is the dominant Ti-species even throughout the LTE 
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atmosphere. 

In the dwarf models, the situation is very different. The higher pressures favor 
molecules over atoms and ions, which impacts the effect of Ti I NLTE on the structure. 
The NLTE T e s = 2700 K model shows increases in the concentrations of TiO + and Ti II for 
r st d < 10 -2 but virtually no change in the concentrations of the dominant TiO and Ti0 2 
molecules. Note the TiO + increases more than the Ti II concentration, indicating that the 
additional Ti II initially created by NLTE overionization of Ti I is converted into TiO + . 

The T eS = 4000 K dwarf model shows a somewhat different behavior. In this model, 
the NLTE effects increase the concentration of Ti II below r st(i ft 10~ 5 significantly so that 
Ti II is the dominant Ti-species below r st( j ft 10~ 6 , whereas in the LTE model TiO becomes 
the dominant species below r st( j ft 10~ 7 . NLTE effects reduce the overall concentration of 
Ti-molecules, but increase the relative importance of TiO + . 

All changes in Ti-molecules occur at small optical depths in the outer regions of the 
model atmospheres. Considering the fact that the molecular lines form deeper in the 
atmosphere than the atomic lines, this indicates that the effects of atomic Ti NLTE on the 
lines of Ti-molecules will be very small. 

3.3. NLTE effects on Ti I line profiles in M Star Spectra 

The influence of NLTE effects on the formation and emergent profiles of the Ti I lines 
can be estimated from the run of the ratio of the line source function Sl to the local 
Planck-function versus the optical depth. This is shown in a series of overview plots in 
Figs. 12 — 16 for all models. The large number of Ti I lines in NLTE masks individual 
transitions, but the figures show the wide range that Sl/ B spans for the different transitions. 
In general, the lines become optically thin for a wide range of optical depths, between 
r std ft 10~ 5 and unity. The electron temperatures and thus the Planck-function drop rapidly 
with smaller optical depths, creating large Sl/B ratios in the outer atmosphere. This 
indicates that Sl has decoupled completely from the thermal pool above the line forming 
region of each line, as expected. 

The Sl/B ratio is overall closest to unity, its LTE value, for the T eff = 3600 K giant 
model. Therefore, this model should show the smallest changes in line profiles of the three 
giant models. However, even in this model many lines have Sl/B ratios very different 
from unity, therefore a number of Ti I lines will be different from the LTE case. For many 
transitions, the ratio Sl/B is less than unity in the line forming region, indicating that the 
absorption lines will show deeper cores than in the LTE case. Note that this effect combines 
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with the effects of the Ti-NLTE on the EOS, which for the models discussed here generally 
reduces the concentration of Ti I, in part canceling the effects due to NLTE scattering. 

We demonstrate the net effect of NLTE on the formation of Ti I lines in giants in 
Figs. 17 to 19. As expected, NLTE effects are smallest for the model with T eS = 3600. In 
the optical spectral region, the changes caused by the Ti I NLTE line formation are very 
small and would be hard to observe due to the enormous crowding of lines in this spectral 
region. NLTE effects for the same Ti I lines are smaller for dwarfs at similar effective 
temperature, cf. Figs. 20 and 21. 

In the cooler models, the Ti I lines form deeper in the atmosphere in a region in which 
the radiation field is nearly Planckian and thus NLTE effects are very small. This is due to 
the enormous background opacity of TiO and other molecules. In the outer atmosphere 
of the cooler models the concentration of Ti I is much smaller than that of TiO, thus the 
effect of the large departures from LTE that we find in these regions on the Ti I line profiles 
is very small. 

In the hotter dwarf models, NLTE effects, in particular on the near-IR lines, are larger. 
The giant models, however, show an opposite behavior (the NLTE effects on the Ti I lines 
are larger for the cooler models). In these models the line forming region of the Ti I lines 
is inside the region in which departures from LTE are significant. In addition, the TiO 
opacity is relatively smaller than in the cooler models. For the Ti I line at A vac « 964lA, 
NLTE effects make the core of the line deeper than the LTE model predicts. This is due to 
line scattering which removes photons from the line core and re-distributes them into the 
line wings. These effects are present in most Ti I lines, but there are a few exceptions in 
the giant model with T e g- = 4000 for which some lines are weaker in NLTE than in LTE (cf. 
Fig. 19). Abundance determinations of Ti, or likely all metals, from near-IR or IR lines 
should therefore include NLTE effects wherever possible. 



4. Summary and Conclusions 

For the model parameters that we have considered so far, effects of Ti I NLTE on the 
TiO bands are very small. This seems to be due to the fact that the line forming region of 
TiO (around r st( j ~ 10 -2 ) is inside of the line forming region of the Ti I lines. Therefore, the 
TiO lines form in a region in which the Ti I atom is basically in LTE and Ti I NLTE effects 
that are important at smaller optical depths cannot affect TiO lines significantly. The 
situation is different for the TiO + molecule which forms from Ti II and O I (cf. Fig. 10). 
This molecule is very sensitive to Ti I NLTE effects and its lines would be helpful indicators 
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of Ti I NLTE effects. 

fn further investigations it will be very interesting to look for NLTE effects in the 
TiO line formation itself. This is certainly feasible with modern numerical techniques once 
adequate data for TiO are available. In addition, NLTE effects of carbon, nitrogen, oxygen, 
Ti I II and, in hotter models, Fe I— II need to be investigated in detail for larger grid of 
models to assess the impact of NLTE on model structures and spectra over wide regions of 
the Hertzsprung-Russell diagram. We are currently preparing such models and will report 
the results elsewhere (Hauschildt et al, in preparation). 
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Fig. 1. — Simplified Grotrian diagram of our Ti I model atom. All 395 levels and 5379 
primary (i.e., full NLTE) lines are shown but the 0.8 million secondary (approximate) NLTE 
lines have been omitted for clarity. 



5. Figures 
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Fig. 2. — Run of the Ti departure coefficients for the T eS = 3200 K, log g = 1.0 giant model. 



Fig. 3. — Run of the Ti departure coefficients for the T eS = 3600 K, log<? = 1.0 giant model. 



Fig. 4. — Run of the Ti departure coefficients for the T e g = 4000 K, log g = 1.0 giant model. 



Fig. 5. — Run of the Ti departure coefficients for the T e s = 2700 K, log g = 5.0 dwarf model. 



Fig. 6. — Run of the Ti departure coefficients for the T eff = 4000 K, log g = 5.0 dwarf model. 
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Fig. 7. — NLTE effects on the Ti ionization and molecule formation for the T eS = 3200 K, 
log g = 1.0 giant model. The thin lines and symbols give the results for the LTE model, 
whereas the thick lines and symbols give the results for corresponding NLTE model. 



Fig. 8. — NLTE effects on the Ti ionization and molecule formation for the T e g = 3600 K, 
logg = 1.0 giant model. The thin lines and symbols give the results for the LTE model, 
whereas the thick lines and symbols give the results for corresponding NLTE model. 



Fig. 9. — NLTE effects on the Ti ionization and molecule formation for the T eS = 4000 K, 
log g = 1.0 giant model. The thin lines and symbols give the results for the LTE model, 
whereas the thick lines and symbols give the results for corresponding NLTE model. 
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Fig. 10. — NLTE effects on the Ti ionization and molecule formation for the T e g = 2700 K, 
log g = 5.0 dwarf model. The thin lines and symbols give the results for the LTE model, 
whereas the thick lines and symbols give the results for corresponding NLTE model. 

Fig. 11. — NLTE effects on the Ti ionization and molecule formation for the T eff = 4000 K, 
logs = 5.0 dwarf model. The thin lines and symbols give the results for the LTE model, 
whereas the thick lines and symbols give the results for corresponding NLTE model. 

Fig. 12. — Line source functions for all Ti I NLTE lines for a giant model with the parameters 
T cfr = 3200K, log<? = 1.0. 

Fig. 13. — Line source functions for all Ti I NLTE lines for a giant model with the parameters 
T cS = 3600 K, logs = 1-0. 

Fig. 14. — Line source functions for all Ti I NLTE lines for a giant model with the parameters 
T cfr = 4000K, logs = 1.0. 

Fig. 15. — Line source functions for all Ti I NLTE lines for a dwarf model with the parameters 
T cfr = 2700K, logs = 5.0. 

Fig. 16. — Line source functions for all Ti I NLTE lines for a dwarf model with the parameters 
T cS = 4000 K, logs = 5.0. 

Fig. 17.— NLTE effects on the Ti I line at A vac « 9641 A and around 5020 A for a giant 
model with the parameters T c s = 3200 K, logs — 1-0- The LTE spectrum (dotted curve) 
uses the same model structure as the NLTE spectrum (full curve) but with all departure 
coefficients set to unity. The fluxes are in arbitrary units. 

Fig. 18. — NLTE effects on the Ti I line at A vac ~ 964lA and around 5020A for a giant 
model with the parameters T c s = 3600 K, logs — 1-0- The LTE spectrum (dotted curve) 
uses the same model structure as the NLTE spectrum (full curve) but with all departure 
coefficients set to unity. The fluxes are in arbitrary units. 

Fig. 19.— NLTE effects on the Ti I line at A vac « 964lA and around 5020A for a giant 
model with the parameters T cS = 4000 K, logs — 1-0- The LTE spectrum (dotted curve) 
uses the same model structure as the NLTE spectrum (full curve) but with all departure 
coefficients set to unity. The fluxes are in arbitrary units. 
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Fig. 20.— NLTE effects on the Ti I line at A vac « 964lA and around 5020A for a dwarf 
model with the parameters T c g = 2700 K, \ogg = 5.0. The LTE spectrum (dotted curve) 
uses the same model structure as the NLTE spectrum (full curve) but with all departure 
coefficients set to unity. The fluxes are in arbitrary units. 



Fig. 21.— NLTE effects on the Ti I line at A vac « 9641 A and around 5020A for a dwarf 
model with the parameters T cS = 4000 K, logg = 5.0. The LTE spectrum (dotted curve) 
uses the same model structure as the NLTE spectrum (full curve) but with all departure 
coefficients set to unity. The fluxes are in arbitrary units. 



